Line Shape Broadening in Surface Diffusion of Interacting Adsorbates 
with Quasielastic He Atom Scattering 
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The experimental line shape broadening observed in adsorbate diffusion on metal surfaces with 
increasing coverage is usually related to the nature of the adsorbate-adsorbate interaction. Here 
we show that this broadening can also be understood in terms of a fully stochastic model just 
considering two noise sources: (i) a Gaussian white noise accounting for the surface friction, and (ii) 
a shot noise replacing the physical adsorbate-adsorbate interaction potential. Furthermore, contrary 
to what could be expected, for relatively weak adsorbate-substrate interactions the opposite effect 
is predicted: line shapes get narrower with increasing coverage. 

PACS numbers: 68.35.Fx, 05.10.Gg, 68.43. Jk 
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The diffusion of atoms, molecules, or small clusters 
on metal surfaces is an elementary dynamical process 
of paramount importance; it constitutes the ground 
step to understanding more complex phenomena in sur- 
face science. Processes such as heterogeneous catalysis, 
molecular-beam epitaxy, crystal growth, chemical vapor 
deposition, or associative desorption are all strongly af- 
fected by the kinetics of diffusion. Different experimen- 
tal techniques have been used to study activated sur- 
face diffusion [H, HI 9 , quasielastic helium-atom scatter- 
ing (QHAS) being one of the most valuable for such a 
purpose 0, H, HJ . In QHAS experiments the observable is 
the so-called differential reflection coefficient. In analogy 
to neutron scattering by liquids Q , this magnitude gives 
the probability that the He-atom beam scattered from 
the diffusing collective reaches a solid angle f2 with an 
energy exchange hw—Ef — Ei and parallel (to the sur- 
face) momentum transfer AK=Ky — Kj, and is propor- 
tional to the dynamic structure factor or scattering law, 
5(AK, to). 5(AK, to) provides information about the dy- 
namics and structure of the adsorbates, thus allowing 
to better understand the nature of adsorbate-substrate 
[B[ and adsorbate-adsorbate Q interactions. In general, 
5(AK,w) usually consists of the quasielastic (Q) peak 
and some weaker peaks attributed to the diffusion process 
and to phonon excitations and adsorbate low frequency 
vibrational modes [namely, the frustrated translational 
(T) modes], respectively. 

Here we present a novel and insightful stochastic ap- 
proach that allows to study and understand in a simple 
manner how the adsorbate coverage influences S*(AK, lu). 
This approach, grounded on the theory of spectral-linc 
collisional broadening developed by Van Vleck and Weis- 
skopf [|| and the elementary kinetic theory of gases Q 
in the Langevin framework, relies on three assumptions: 
(i) the adsorbate-substrate interaction is described by 
a (deterministic) adiabatic adsorption potential, (ii) the 



temperature dependent nonadiabatic coupling to the sub- 
strate electronic and vibrational excitations is accounted 
for by a Gaussian white noise, and (iii) adsorbate- 
adsorbate interactions are described by a shot noise 
[Icl Hl| , as proposed in Ref. UH. We call this approach 
the interacting single adsorbate approximation. Gaussian 
white noise has been used to characterize the rich vari- 
ety of different transport regimes in surface diffusion [l3| . 
Within a similar context, the shot noise has been used 
to study thermal ratchets 11411 . m ean first passage times 
[l5j |. and jump distributions . 

To explain the line-shape broadening undergone by 
the Q peak with increasing coverage (#), high demand- 
ing molecular dynamics calculations within the Langevin 
framework have been carried out Though the cal- 

culations could reproduce the trend with 0, the dipole- 
dipole interaction potential considered was not able to 
correctly provide the experimental broadening observed. 
As we show here, this (9-dependent behavior can also be 
reproduced by a simple collision model, with collisions 
being described by a shot noise. Note that, providing 
collective behaviors are not relevant (as happens, for in- 
stance, at low temperatures jl8|). adsorbates undergo a 
relatively large number of collisions for relatively long 
time scales reaching the statistical or stochastic limit. 
This proves that the nature of the adsorbate-adsorbate 
interaction could not be very relevant regarding the line 
shape broadening (as well as its trend). As seen here, the 
broadening is strikingly related to the surface corrugation 
(which strongly couples diffusion and low frequency vi- 
brations) as well as to the friction associated to increas- 
ing 9. This important result emerges when combining 
the numerical simulations using our model with a fully 
analytical treatment, which allows to understand and to 
interpret the process in terms of two limiting cases: run- 
ning or diffusive and bound trajectories [l9| . 

The starting point of our approach consists of express- 
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ing the dynamic structure factor as 0, HO] 

S(AK,w) = J e-™* I(AK,t)dt, (1) 

where 

J(AK,t) = ( e - iAK -[ R (*)- R (°)]) = ( e -i&Kftv^ K (t')dt'^ 

(2) 

is the intermediate scattering function. In the brack- 
ets denote an ensemble average and «ak is the adparticle 
velocity projected onto the AK direction (AK = || AK||). 
After a second order cumulant expansion in AK in the 
right-hand side (rhs) of the second equality of ([2]), one 
obtains 

/(AK,f)«e- Aif2 /o(t-*')C(t')dt' ) ( 3 ) 

where C(t) — (wak(O)wak(O) 1S the velocity autocorre- 
lation function. This is the so-called Gaussian approxi- 
mation Q, which is exact when velocity correlations of 
third or higher order are negligible. Despite its limita- 
tions, it provides much insight into the dynamical process 
by allowing an analytical treatment of the problem. 

Determining the line shape of S(AK,ui) through ^ 
thus requires us to simulate the adsorbate dynamics. 
This is done within a purely Langevin framework that in- 
cludes the three elements mentioned above. Considering 
one dimension for simplicity, the motion of an adsorbate 
interacting with another adsorbates and a temperature- 
dependent substrate is ruled by the generalized Langevin 
equation [11] 

x(t) = - [ r)(t-t')x(t')dt' + F[x(t)]+SR(t), (4) 
Jo 

where x is the adsorbate coordinate, rj(t) is the memory 
function, F——W is the deterministic force per mass 
unit (V is the periodic adsorbate-surface interaction po- 
tential, with period a), and SR(t)=SRc(i)+SRs(t) is the 
noise source (G and R refer to Gaussian white noise and 
shot noise, respectively). 

Gaussian white noise is defined by (Ra{t))—0 and 
{R G (t)R G (t'))=2m~fk B T5(t' - t), where m is the adsor- 
bate mass, T is the surface temperature, and 7 is the 
(constant) friction coefficient measuring the strength of 
the adsorbate-substrate coupling. On the other hand, the 
shot noise term in our model is Rs(t)=^2k=i bk{t — tk), 
where bk(t — tk)=Ck\'e~ x ( t_ * fc ) is the pulse or impact 
shape, with t — tk>0 and Ck giving the intensity of the 
collision impact. The pulse shape indicates that col- 
lisions are assumed to be sudden (strong but elastic) 
and after-collision effects relax exponentially at a con- 
stant rate A'. The probability for K collisions to hap- 
pen after a time T follows a Poisson distribution [1 01 ] . 
P K (T)=(K\y 1 (XT) K er XT , where A means the average 
number of collisions per time unit. Since collisions take 
place randomly at different orientations and energies, 



we reasonably assume g(ck)=a x e Cfc / Q , with Cfc>0 and 
a=sJm/k B T [H. 

The rate A' defines a decay time scale for each colli- 
sion event, r c =l/A'. If t c is relatively small (collision 
effects relax relatively fast), the memory function associ- 
ated to the shot noise in ((4]) is also local in time. Then, 
r](t — t'}~n 6(t — t') and, resorting to the Markovian ap- 
proximation [10J, Eq. ((4| reduces to 

x(t) = -rjx(t) + F[x(t)} + SR(t), (5) 

with a total constant friction n=j+\ and where 
(6R(t)6R(t'))=2mr)k B T5(t' - t). To estimate the value 
of A, the elementary kinetic theory of transport in gases 
Q can be considered, where diffusion is proportional to 
the mean free path I, and the latter is proportionally 
inverse to both the density of gas particles and the effec- 
tive collision area when a hard-sphere model is assumed. 
By means of simple arguments [12j , it is possible to find 
an also simple relationship between A and 9, given by 
\=(&p9 / a 2 )yJk B T /m. Accordingly, increasing 9 (and/or 
T) means an also increase of the collision frequency. 

In Fig. [1] we have plotted 2D calculations for S(AK, us) 
with two different types of corrugation: (a) low corru- 
gation, assuming V=0, and (b) a non separable corru- 
gated potential [5(]. We have considered Na on Cu(001) 
at A!^=1.23 A -1 , surface temperature of 200 K and 
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FIG. 1: (Color online.) Q peak for Na on Cu(001) at 
AA"=1.23 A -1 , T=200 K, and 01=0.028 (blue dashed line) 
and 02=0.18 (red solid line), (a) Flat surface potential (V=0) 
and (b) non separable corrugated potential [5|. 
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two coverages, #i=0.028 and #2=0. 18, as in the exper- 
iment 0. After calculating 2D trajectories from Eq. ([5]), 
7(AK, t) is obtained using the rhs of the first equality in 
Eq. © and fitted to a certain analytical function (see be- 
low). Finally, the Fourier transform of that function gives 
5(AK, u). As is clearly seen in Fig.QJa), the Q peak gets 
narrower as 9 increases for low corrugation. On the con- 
trary, the presence of a potential leads to the broadening 
observed in Fig.QJb). For instance, the full widths at half 
maximum (FWHM) for three coverages at AA"=1.26 A -1 
are: r=120, 190, and 240 ^eV for (9=0.028, 0.106, and 
0.18, which are in good agreement with the experimen- 
tal ones, r=110, 150, and 270 ^eV, respectively; these 
widths differ from those found in Ref. [6j (r=110, 130, 
and 160 ^eV, respectively) by using a dipole-dipole in- 
teraction in Langevin molecular dynamics simulations. 

To understand our results we can resort to two (analyt- 
ical) limit cases: pure diffusion and anharmonic oscilla- 
tions. After assuming V—Q in Eq. ([5|), one easily obtains 
C{t)^(vl)e-^ (where (v%)=k B T/m), which leads to 



J(AK, t) = cxp [-x 2 (e~* + nt- 1)] 



(6) 



after using Eq. © with x 2 ={vl)^K 2 /n 2 . On the 
other hand, for a damped harmonic oscillator (l9l . 
[H we have C(t )=(^) e - ? ' t / 2 cos{Qt + <5)/cos<5, with 
lu=^/lu 2 — (ri/2) 2 , <5=(tan) -1 (?7/2(D), and uj being the 
harmonic frequency. Again, using Eq. ([3]), we obtain 



I(AK, t) = exp 



2 1 

XT 

LUUJo 



cos S - e 1,4/2 cos(u;t - S) 



(7) 

Unlike (J6]), Eq. |7]) displays an oscillatory but exponen- 
tially damped behavior which does not decay to zero. 
Nevertheless, it approaches © in the limit lo 0. Here 
we are not interested in the damped harmonic oscilla- 
tor, but in an anharmonic one, which arises when the 
parameters in ||7J) are left free. 

If deviations from the Gaussian approximation are not 
very important, 7(AK, t) will display features typical of 
the behaviors described by both ([6]) and ([7|). Indeed, 
even if the Gaussian approximation does not hold, one 
can still make use of a working formula in order to ex- 
tract relevant information about the process from the 
numerical results. Trajectories display both temporary 
trapping in potential wells and periods of time where the 
flight is unbound. Therefore, assuming a model of run- 
ning (i?) and bound (B) trajectories, one can consider 
C{t)=aCn{t) + (l — a)Cs{t), where Cr and Cb correspond 
to the velocity autocorrelation functions for a flat sur- 
face and a damped anharmonic oscillator, respectively. 
7(AK,i) can be then written as 



/(AK,t) « [I R (AK,t)] a [I B (AK,t)} 



(8) 



where Ir and 1b are given by Eqs. © and 0, respec- 
tively. This expression, with free parameters, is used as 



our working formula and allows for a distinction between 
the contributions arising from the unbound or diffusive 
motion and the bound or vibrational one. As seen in 
Fig. [21 Eq. ([5]) fits fairly well the Langevin numerical re- 
sults obtained for the corrugated surface potential [see 
Fig. [ljb)] and both coverages. The fitting parameters 
compared with the nominal values are given in Table [H 
though fitted and nominal values are different, their or- 
der of magnitude and trend are correct, regarding such 
differences to deviations from the Gaussian approxima- 
tion. More importantly is, however, the resulting value 
for a: a x =0.0294 for 6 X , and a 2 =0.0410 for 2 . From 
these values, it is clear that the contribution to I(AK, t) 
primarily comes from bound motions. Note from Q that 
running trajectories lead to a relative much faster decay 
of J(AK, t) than the bound motion [see Eq. (J7J)], which 
delays such a decay. This is so because the bound mo- 
tion keeps correlations for times longer than the diffusive 
one, since the latter provokes a fast dephasing of the (cor- 
relation) oscillating terms that appear in the rhs of the 
first equality in Eq. |2]). In this way, this explains, first, 
that S(AK,w) is about two orders of magnitude broader 
in the flat case than in the corrugated one (see Fig. [l| . 
And, second, since a 2 >ai there is a slightly larger frac- 
tion of running trajectories for 02, which leads to an also 
slightly faster decay of 7(AK,f), and therefore to ob- 
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FIG. 2: (Color online.) 7(AK,t) for Na on Cu(001) at AK = 
1.23 A" 1 , T=200 K, and (a) <9i=0.028 and (b) 2 =O.18. Open 
circles indicate the numerical values obtained from the simu- 
lation and solid lines are the numerical fitting using Eq. (JHJ). 
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coverage 


uj (a.u.) 


r\ (a.u.) 
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X 


^nom 

ef 


2.20xl0 -4 
1.39xl0" 4 


2.54x10" 
3.52x10" 


5 
-5 


0.058 
0.184 


3.15 
4.64 


nnom 

ef 


2.19xl0" 4 
1.37xl0 -4 


4.69x10" 
4.28x10" 


6 
-5 


0.107 
0.276 


1.71 
2.82 



TABLE I: Nominal and fitted values involved in the calcula- 
tion of I(AK,t) for 0i=O.O28 and 2 =O.18. 



serve broadening in S(AK,w) with increasing 9. It is 
also worth stressing that, with increasing coverage, dif- 
fusion decreases (according to Einstein's relation, it goes 
like f7 _1 ) and the running trajectories display short flights 
in small surface regions. 

To conclude, the simple stochastic model together with 
the theoretical analysis presented in this work shows that 
the experimental broadening observed in the Q peak 
with increasing 9 (from low to intermediate regimes) 
is related to the bound motion undergone by the ad- 
sorbates inside the potential wells, and not only to 
purely diffusive motion or to the particular form assumed 
for adsorbate-adsorbate interactions. Moreover, regard- 
less the adsorbate-substrate interaction potential, in our 
model the line shape of S(AK,w) is also related to rj, 
which gathers the thermal effects caused by the surface 
(7) and the collisions among adsorbates (A). Since 7 
is assumed to be constant, the broadening is thus di- 
rectly related to an increasing friction A as 9 increases. 
The idea of replacing the dipole-dipole interaction by a 
shot noise is crucial because, for long time processes with 
a high number of collisions, the statistical limit seems 
to wipe out any trace of the true interaction potential. 
As happens with the models proposed in the literature 
[S H EH' ours gi yes a smaller diffusion than expected 
with increasing 9, which could be related to the fact that 
9 might also influence the type of adsorbate-substrate in- 
teraction, and therefore modify parameters such as the 
activation barrier or the surface friction in our model. 
Thus, in Ref.Qil, the Lau-Kohn long range interaction via 
intrinsic surface states is proposed to use a ^-dependent 
surface electronic density modifying the corrugated po- 
tential and to explain the remarkable increasing of the T 
mode frequency with coverage at 50 K which a dipole- 
dipole or substrate mediated coupling cannot explain it. 
This increasing is also observed at 200 K and parallel 
momentum transfer of —2 A" 1 to be 11%, which is sim- 
ilar to our value (around 7%) when going from 0=0.028 
to 0.106 at 2.10 A -1 . Although further investigation at 
microscopic level and calculations from first principles 



are needed, our simple stochastic model at low coverages 
and surface temperatures is also able to provide a comple- 
mentary view of diffusion and low-frequency vibrational 
motions, described by the peaks around or near zero en- 
ergy transfers in the scattering law. 
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